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Abstract. Studies of the diffuse X-ray emitting gas in galaxy clusters have provided powerful con- 
straints on cosmological parameters and insights into plasma astrophysics. However, measurements 
of the faint cluster outskirts have become possible only over the last few years. Here, we present 
results from Suzaku observations of the Perseus Cluster, which provide our best measurements of 
the thermodynamic properties of the ICM at large radii to date. In particular, we focus on the details 
of the data analysis procedure and discuss the evidence for a clumpy distribution of the gas in the 
outskirts, which is important for understanding the physics of the ongoing growth of clusters from 
the surrounding cosmic web. 
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MOTIVATION 

In the hierarchical picture of large-scale structure formation in the Universe, clusters 
of galaxies, which are the most massive objects, are also the latest to form. Galaxy 
clusters are enormous knots in the cosmic web, located at the intersections of large- 
scale structure filaments, from which they are still accreting matter at the present time. 
Detailed observations of cluster outskirts therefore allow us not only to determine the 
cluster properties accurately but also to witness and understand large-scale structure 
formation as it happens. 

The positions of clusters at the high end of the mass spectrum makes them particularly 
useful and sensitive cosmological probes, provided that their total mass can be measured 
accurately [for a review, see 1]. Hydrostatic mass measurements extending to large radii, 
used in combination with other thermodynamical X-ray observables and measurements 
of weak-lensing, galaxy velocity dispersion, and the Sunyaev-Zel'dovich (SZ) effect, 
provide critical constraints on structure formation models and a robust comparison point 
for simulations. 

Observations of cluster outskirts can be used to search for the boundary between 
virialized and infalling material, beyond which the hydrostatic equilibrium assumption 
breaks down. At larger radii, we can therefore expect to detect inhomogeneities in the 
ICM, such as infalling matter clumps and signatures of accretion shocks, which could 
provide a direct view of the physical processes by which clusters grow. Given even 



a limited resolution in the azimuthal direction we can furthermore ask whether the 
signatures of gas accretion are ubiquitous or occur preferentially along the major axis or 
in the direction of known large-scale structure filaments in supercluster environments. 

THE ROLE OF THE SUZAKU SATELLITE 

Until recently, only a limited amount of information regarding the X-ray surface bright- 
ness of cluster outskirts was available, primarily based on ROSAT observations [2, 3]. 
The low surface brightness signal and relatively high instrumental noise of the detectors 
onboard Chandra and XMM-Newton have limited robust measurements of the tempera- 
ture of the hot, diffuse intracluster medium (ICM) to radii less than approximately half 
of the virial radius [e.g. 4]. Implicitly, measurements of the total cluster mass, deter- 
mined under the assumption of hydrostatic equilibrium, could not be extended outside 
half of the virial radius except with the use of parametrized extrapolations. 

The lower and more stable particle background of Suzaku has enabled a breakthrough 
in determining the ICM properties out to the virial radius in several bright, relatively 
nearby clusters [5, 6, 7, 8, 9, 10]. In all cases, the Suzaku results point towards temper- 
ature profiles which decline smoothly from a few hundred kpc out to large radii, which 
is qualitatively in agreement with predictions from numerical simulations [e.g. 11]. For 
several clusters, it was found that, at large radii, the entropy profile flattens compared 
to the power-law increase with radius observed in the central regions [5, 10], possibly 
pointing towards the onset of non-hydrostatic equilibrium effects and convective insta- 
bility around and beyond the virial radius. 

THE OUTSKIRTS OF THE PERSEUS CLUSTER: AN IDEAL 

TARGET FOR SUZAKU 

The Perseus Cluster of galaxies (A426) is the brightest, extended X-ray source. It is both 
closer (at z=0.0183) than any of the clusters for which Suzaku has been previously used 
to measure the ICM properties at large radii, and it has a significantly higher X-ray flux. 
The large angular size of the cluster mitigates the need to precisely model the effects of 
Suzaku's complex point spread function (PSF), making the Perseus cluster an ideal target 
to study cluster outskirts. Furthermore, the large angular size permits the use of larger 
extraction regions for the same fractional interval of the virial radius, which translates 
into collecting signal faster and reducing the expected root mean squared variations of 
the background point source flux. A wealth of information is already available about 
the central parts of the Perseus cluster, from the Uhuru and Einstein satellites up to a 
recent 1.4 megasecond Chandra image [12]. This information, in combination with the 
unprecedented Suzaku measurements of the temperature out to the virial radius which 
will be presented in the following sections, provides the most detailed view of a galaxy 
cluster ever achieved. 



FIGURE 1. The X-ray surface brightness map of the NW and E arms of the Perseus Cluster, superim- 
posed to scale on an optical SDSS image. The dashed circle marks r200. The inset shows a deep Chandra 
image of the cluster's bright central region. 

The pilot project on the outskirts of the Perseus Cluster 

Between 2009 July 29 and 2009 August 22, Suzaku observed a mosaic of 14 fields 
placed along two arms towards the east and northwest from NGC1275, the central galaxy 
of the Perseus cluster. NGC1275 is monitored regularly as a calibration source, therefore 
in addition we used one central pointing performed on 2009 August 26, which is as close 
in time as possible to our mosaic fields and was available from the public archive. We 
focus on the data obtained with the three available X-ray imaging spectrometer (XIS) 
cameras. The main results on this study are described in [13]. The purpose of the current 
manuscript is to complement the main journal article and supply additional information 
regarding the details of the data analysis performed. 

The data were reduced using the tools available in the HEAsoft package (version 
6.9) to create a set of cleaned event lists with hot or flickering pixels removed. All the 
standard recommended screening criteria were applied (see e.g. http : / /heasarc . 
nasa.gov/docs/suzaku/processing/criteria_xis.html) and we se- 
lected times with the geomagnetic cut-off rigidity COR > 6 GV. We also verified using 
the method described in [14] that the geometry of Earth's magnetic field was not favor- 
able for generating strong solar wind charge exchange (SWCX) signatures during our 
observations. 

The X-ray surface brightness image, extracted in the 0.7-7 keV energy band and 



corrected for vignetting and instrumental background, is shown in two false-color strips 
in Fig. 1. We visually identified bright point sources marked by red circles in the figure. 

Background analysis 

The non X-ray background (NXB) spectrum was constructed from the dark Earth 
database using the standard method in which the COR distributions of the on-source 
and the background data are matched. The cosmic X-ray background (CXB) is typically 
modeled with three components: a thermal component with kT~ 0.09 keV to account 
for the local hot bubble (LHB) emission, a second thermal component with kT~ 0.2 keV 
for the Galactic halo (GH), and a power-law with T = 1.41 to account for the integrated 
emission of unresolved point sources [e.g. 15, 16]. Additionally, we include a thermal 
component with kT~ 0.6 — 0.8 keV (hot foreground, HF), which is sometimes required 
to describe the Galactic foreground, particularly at lower Galactic latitudes. 

Data from the Suzaku mosaic beyond the estimated virial radius of Perseus (the outer 
two pointings along each arm), as well as from the ROSAT All Sky Survey (RASS) can 
be used to constrain the parameters of the background model and estimate their spatial 
variations. From RASS, we used 5 circles with a radius of 1 degree located 3° 10' to 
the N, NE, E, SE, and S of the cluster center, avoiding contamination from Algol and 
AWM7 on the western side. Solar abundances were assumed for all thermal components 
[17]. Table 1 shows an excellent agreement between ROSAT and Suzaku for all the 
model parameters. For our final CXB model, we chose the best- fit parameters of the 
GH, HF, and power-law based on the Suzaku data. The LHB parameters, which cannot 
be constrained with Suzaku, were fixed based on the ROSAT data. 

To determine the spatial variations of the Galactic foreground, we fitted spectra from 
the RASS circles leaving the normalizations of the GH, HF and LHB components to vary 
independently for all 5 spectra. The temperatures of all thermal components were fixed 
to our choice of the final CXB model parameters described above. The LHB contributes 
very little to the emission above 0.5 keV, thus we do not calculate the effect of systematic 
uncertainties on the LHB parameters. To estimate the systematic uncertainties on the 
CXB power-law, we fitted 8 Suzaku background spectra (each of the outer two pointings 
along the two arms was divided in half) individually in the 2-7 keV band and computed 
the variance of the 8 values of the normalization. The area of each of the 8 regions used 
was approximately 160 sq arcmin, thus for smaller spectral extraction regions of area A 
the variance is expected to increase by a/ 160 arcmin 2 /A. 



Stray-light rejection for the spectral analysis 

For each annulus, we calculated the difference between the observed spectrum and 
the assumed CXB model presented above. This yields a background-subtracted cluster 
signal. We then divided by this signal the corresponding stray light spectrum extracted 
from the event files generated by our ray-tracing simulations. We included in the spectral 
fitting only the parts of each observed spectrum between 0.7-7 keV where the contami- 



TABLE 1. Cosmic X-ray Background parameters. The temperatures are given in keV, 
all normalizations represent the standard XSPEC value normalized to 10~ 3 per 20 2 n sq 
arcmin area. Errors are given at the AC=1 level. 



Parameter Suzaku ROSAT Systematic uncertainty 



LHB 


kT 


0.091 (fix) 


0.091 ±0.003 






norm 


1.70 (fix) 


1.70 ±0.07 




GH 


kT 


0.18 ±0.02 
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HF 
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0.50 ±0.07 


0.31(S)-0.60(SE) 


Power-law 
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1.41 (fix) 


1.41 (fix) 






norm 


1.096 ±0.016 


1.138 (fix) 


9.4% xy/ 160 arcmin 2 /A 
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FIGURE 2. Effects of systematic uncertainties in the CXB parameters on the cluster normalization and 
projected temperature and metallicity profiles. 

nation was below 20%. Because the stray light spectrum peaks between 1-1.2 keV, even 
in regions with high stray-light contamination the spectral contribution of the stray light 
above 1.5 keV becomes typically less than 10%, and without exception less than 20%. 

Results 

Combining our results with those obtained from an ultra-deep Chandra observation of 
the cluster center[18], we find that the Suzaku and Chandra radial profiles show excel- 
lent agreement where they intersect, and together measure the temperature and metal- 
licity structure of the intra-cluster gas with high precision and previously unattained 
spatial resolution out to the virial radius (/200) [13]. In the narrow interval spanning 



0.95-1.05r2oo> the temperature is approximately a third of the peak temperature. The 
systematic errors introduced by the CXB modeling are typically smaller than the statis- 
tical la errors, as demonstrated in Fig. 2. 

We are also able to obtain exquisite deprojected radial profiles of the electron density, 
entropy, and pressure. Assuming that the total mass distribution follows an NFW profile, 
as suggested by numerical simulations [19], we used the data from the dynamically 
relaxed northwestern arm to determine the best-fit total mass profile and obtain an 
unprecedently accurate gas mass fraction profile, which extends out to the virial radius 
of the Perseus Cluster (shown in the right panel of Fig. 3). 

This profile reveals no evidence for the puzzling deficit of baryons at r > 0.5r2oo 
inferred from some previous studies of other systems, using lower quality data and/or 
extrapolated models [e.g. 20, 4]. At r 50 o, the cluster has the expected (approximately 
Universal) baryon fraction. Within r < 0.5r2oo> then, the physics of the X-ray emitting 
gas appears relatively simple and X-ray measurements can be used robustly for cosmo- 
logical work. At larger radii, the apparent gas mass fraction exceeds the cosmic mean 
baryon fraction measured from the CMB [21]. This excess also correlates with an ap- 
parent flattening of the entropy profile. These results are robust against all known sys- 
tematic uncertainties. The most plausible explanation for these properties is that the gas 
at large radii is clumpy 1 . The deprojected radial profile of the electron density shows ex- 
cellent agreement with previous ROSAT measurements extending out to ~ 1.4 Mpc[22], 
therefore we can exclude any residual stray light contamination in the Suzaku data as a 
possible reason for the boosted electron density at large radii. 

Outside the central region, and inside the radius where clumping becomes important, 
the measured f gas profile shows good agreement with recent numerical simulations [25], 
where a semi-analytic model was used to calculate the energy transferred to the intra- 
cluster gas by supernovae and AGN during the galaxy formation process. Extrapolating 
this model into the outskirts where clumping is important, we used its predictions to- 
gether with the measured f gas ( r ) to determine by how much the electron density must 
be overestimated to produce the difference between the data and the model. Correcting 
the electron density using this factor, we find that the entropy profile becomes consistent 
with the power-law profile expected from models of gravitational structure formation 
[23]. Additionally, the pressure becomes consistent with that expected by extrapolating 
the average profile of a sample of clusters previously studied with XMM-Newton[24]; 
this "universal" pressure profile is also motivated by numerical simulations [26]. Cor- 
recting for the bias introduced by gas clumping is the only straightforward way to bring 
both the entropy and pressure profiles in agreement with the expected models. 

Numerical simulations predict gas clumping in the cluster outskirts [27, 28], but the 
amount of clumping depends on many physical processes in the ICM which are currently 
uncertain (e.g. viscosity, conduction, star formation feedback, magnetic fields), and on 
the numerical schemes used [29]. To explore the implications of gas clumping for cluster 
astrophysics, we must combine efforts in both observation and theory. 



if the density is non-uniform, then the average of the square of the electron density, which determines the 
X-ray bremsstrahlung emission, can be larger than the square of the average electron density, (n%) > (n e ) 2 . 
Determining the electron density from X-ray observations as \f (nf) will overestimate the true n e . 




FIGURE 3. Left: The deprojected electron density (n e ), entropy (K), and pressure (P) profiles measured 
with Suzaku along the NW arm [13] are shown with red data points. The top panel shows the agreement 
with the n e measured with ROSAT in the N quadrant [22]. The red line shows the same profiles corrected 
for clumping. The expected entropy profile from simulations of gravitational collapse[23] is a power-law 
with index j3 ~ 1.1, over-plotted as a black dotted line in the entropy panel. The average profile of a 
sample of clusters previously studied with the XMM-Newton satellite within ^0.5/"2oo[24] is shown with 
a solid black curve in the pressure panel; its extrapolation to r?oo is shown with a dotted black line. Right: 
The integrated, enclosed gas mass fraction profile for the NW arm. The cosmic baryon fraction from 
WMAP7[21] is indicated by the horizontal solid black line; accounting for 12% of the baryons being in 
stars gives the expected fraction of baryons in the hot gas phase, shown as a dashed black line. Predictions 
from numerical simulations [25] are shown in blue. The bottom panel shows the inferred clumping factor. 



THE PERSEUS CLUSTER KEY PROJECT AND FUTURE 

OUTLOOK 

Numerical simulations predict large variations in the gas clumping factor as a function 
of azimuth and dynamical state of the ICM. It is therefore crucial to have measurements 
of similar precision to our pilot project along other directions in Perseus, as well as in 
other systems. The Perseus Cluster Key Project will extend the existing mosaic by 6 
additional arms (2 in AO-5 and 4 in AO-6), providing a uniform sampling in azimuth. 
In addition, Key and Large projects covering the Coma Cluster (a very massive, but 
unrelaxed system) and A2199 (a less massive, relaxed cluster) will provide crucial 
information about the variation of the thermodynamic and f gas profiles and of the 
clumping factor in systems of different masses and dynamical states. High-resolution 
imaging with Chandra will further reveal the shapes and size distribution of the brightest 
gas clumps. 
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